8 Fisher, Scharf, & Lahav 



Acknowledgments 

We thank A.J.S. Hamilton, P. Lilje, D. Lynden-Bell and H. 
Mo for helpful comments and discussions. We also thank M. 
Strauss, J. P. Huchra, M. Davis, and A. Yahil for permission 
to use the 1.2 Jy IRAS survey data prior to its publication. 
KBF and CAS acknowledge a SERC postdoctoral fellowship 
and studentship respectively. 



REFERENCES 

Bertschinger, E. 1992, in New Insights into the Universe, ed. V. 

J. Martinez, M. Portilla & D. Saez (New York: 

Springer- Verlag), in press 
Bertschinger, E. & Dekel, A. 1989, ApJL, 336, L5 
Dekel, A., Bertschinger, E., & Faber, S.M. 1990, ApJ, 364, 349 
Dekel, A., Bertschinger, E., Yahil, A., Strauss, M.A., Davis, M., 

&z Huchra, J. P. 1993, ApJ, in press 
Efstathiou, G., Bond, P., & White, S.D.M. 1992, MNRAS, 258, 

IP 

Feldmann, H., Kaiser, N., &z Peacock, J. 1993, preprint 
Fisher, K. B. 1992, Ph.D. thesis. University of California, 
Berkeley 

P^isher, K. B., Davis, M., Strauss, M. A., Yahil, A., & Huchra, J. 
1^ P. 1993a, ApJ, 402, 42 

Fisher, K. B., Davis, M., Strauss, M. A., Yahil, A., & Huchra, J. 
QQ P. 19936, MNRAS, in press 

^Frenk, C. S., White, S. D. M., Efstathiou, G., & Davis, M. 1990, 
ApJ, 351, 10 

Gorski, K., Davis, M., Strauss, M. A., White, S. D. M., & Yahil, 
t ■ A. 1989, ApJ, 344, 1 
CN Hamihon A. J. S. 1992, ApJL, 385, L5 
OHamihon A. J. S. 1993, ApJL, 406, L47 

On Hamihon A. J. S., Kumar, P., Lu, E., & Matthews, A. 1991, 

O ApJL, 374, LI 

CO Kaiser, N. 1987, MNRAS, 227, 1 

On Kaiser, N., Efstathiou, G., Elhs, P., Frenk, C, Lawrence, A., 
^ Rowan-Robinson, M., & Saunders, W. 1991, MNRAS, 252, 1 
^ K aiser, N. &z Lahav, O. 1988, in Large-Scale Motions in the 

I Universe, ed. V. C. Rubin, G. V. Coyne, S. J. (Princeton: 
^ Princeton University Press), p. 366 

Lahav, O., Yamada, T., Scharf, C, &l Kraan-Korteweg, R. C. 
1993, MNRAS, 262, 711 
^ Lilje, P. & Efstathiou, G. 1989, MNRAS, 236, 851 

Lynden-BeU, D., Lahav, O., & Burstein, D. 1989, MNRAS, 234, 
677 

McGiU, C. 1990, MNRAS, 242, 428 

Regos, E. & Szalay, A. S. 1989, ApJ, 345, 627 

Rowan-Robinson, M. et al. 1990, MNRAS, 247, 1 

Peebles, P.J.E. 1980, The Large-Scale Structure of the Universe, 

(Princeton: Princeton University Press) 
Peebles, P.J.E. 1973, ApJ, 185, 413 

Saunders, W., Rowan-Robinson, M., Lawrence, A., Efstathiou, 
G., Kaiser, N., Ellis, R.S., & Frenk, C.S. 1990, MNRAS, 242, 
318 

Scharf, C, Hoffman, Y., Lahav, O., & Lynden-BeU, D. 1992, 

MNRAS, 256, 229 
Scharf, C. & Lahav, O. 1993, MNRAS, in press 
Strauss, M.A., Davis, M., Yahil, A., & Huchra, J.P. 1990, ApJ, 

361, 49 

Strauss, M. A., Huchra, J. P., Davis, M., Yahil, A., Fisher, K. B. 

& Tonry, J. 19926, ApJS, 83, 29 
Strauss, M. A., Yahil, A., Davis, M., Huchra, J. P., & Fisher, 

K. B. 1992a, ApJ, 397, 395 
Yahil, A., Strauss, M.A., Davis, M., & Huchra, J.P. 1991, ApJ, 

372, 380 



Spherical Harmonics and Redshift Distortion 7 



et al. (1991) found /? = 0.91q'j5 (I''') in ^ comparison of 
the predicted IRAS peculiar velocity field with direct mea- 
surements using the deeper, albeit sparser, QDOT redshift 
survey. Measurements based on the predicted acceleration 
of the Local Group using IRAS catalogs have also yielded 
high estimates 0.4 < /? < 1.0 (Rowan- Robinson et al. 1990; 
Strauss et al. 1992a). Hamilton's (1993) analysis of redshift 
distortions in the correlation function for the 1.936 Jy IRAS 
survey gave a slightly lower value, /? = 0.661q'22 
though consistent with the SHA to within the quoted la 
error. 

These results provide strong evidence of a high, near clo- 
sure, value of Qo on large scales the bias factor for IRAS 
is near unity. However, one should be cautious in claiming 
that these results provide definitive proof of Qo = 1. Sim- 
liar estimates of /? using optical surveys have yielded lower 
values (0.25 < /? < 0.65) for the optical dipole (Lynden- 
Bell, Lahav, & Burstein 1989) and (/? ~ 0.5) in comparisons 
of optical samples with peculiar velocities (Hudson 1993, in 
preparation)) suggesting that at least the ratio of the IRAS 
and optical bias factors may differ from unity. 

The estimate of /? obtained with the SHA is largely free 
of systematic errors. The main systematic error arises from 
uncertainties in the shape of the IRAS power spectrum on 
large scales; analyses of the extended QDOT redshift sur- 
vey (Saunders et al. in prepartion) should greatly reduce 
this uncertainty. The SHA method has the further advan- 
tage of requiring only redshift data. This should be con- 
strasted with the POTENT technique which, although con- 
ceptually very elegant, can be affected by systematic errors 
in the direct measurements of peculiar velocities (e.g., the 
Malmquist bias). Unlike dipole analyses of the Local Group 
acceleration which are plagued by systematic errors arising 
from redshift distortions (the "rocket" effect discussed in 
§ 2), the SHA treats the redshift distortion in an internally 
self-consistent way. Estimates of /? from the redshift cor- 
relation function (Hamilton 1993) are prone to system- 
atic errors arising from non-linear distortions on small scales 
and are extremely noisy on larger scales where the signal to 
noise ratio in the estimate of the correlation function is low. 
Morever, linear analyses based on Kaiser's formulation of 
the distortion neglect the term (2 -|- — Xai^j ^ 

(Kaiser 1987, equation 3.3); although the relative impor- 
tance of this term falls with distance, it should be included 
when analysing the redshift space clustering of local, and 
usually well sampled, structure. 

We are currently examing the question of how to de- 
sign an optimal weighting scheme that will yield a minimum 
variance estimate of /? using techniques similar to those dis- 
cussed in Strauss et al. (1992a). Several main criteria should 
be considered when picking the functional form of the differ- 
ent weighting functions. Firstly, the derivative of the win- 
dow function should be nonvanishing since the distortion 
indentically vanishes when df / dr = 0. Secondly, the various 
weighting functions should sample different regions of space 
to avoid a high covariance between the structure sampled. 
Thirdly, the windows defined by the weighting functions 
should not give excessive weight to distant galaxies which 
will be poorly sampled in current redshift surveys. Fourthly, 
the windows should sample a large enough volume to ensure 
that the harmonics are indicative of their rms values. One 



should keep in mind, however, that the statistical error for 
the weighting scheme used in our analysis was comparable 
to the systematic error induced by uncertainties in the shape 
of the power spectrum. Consequently, further optimization 
will improve the quoted errors only if coupled with more 
reliable estimates of the power spectrum on large scales. 

The SHA formalism for the redshift distortion in galaxy 
samples is also very convenient for comparison with two 
other major probes of the cosmic structure, peculiar veloci- 
ties and the microwave background. The difference between 
the real and redshift space harmonics can be written (cf.. 
Equations (9), (9), and (11)) for I > 2 as 



fii ' - / dk PM(k) 



dr 4>{r)'^^^ j'i{kr) 



(18) 



where PM{k) = PR{k)/b^ is the power spectrum of the mat- 
ter fluctuations. 

This is entirely analogous to expressions for bulk flows 
and peculiar velocities. In particular, it is similar to the ex- 
pansion of the radial velocity field in harmonics (Regos & 
Szalay (1989), equations 44-45) : 



dk PM{k) 



d'^r ^{r)j'i{kr) 



(19) 



where $(s) is the generalized (anisotropic) selection function 
for the observed radial velocities. A comparison of Equa- 
tions (18) and (19) shows that the redshift distortion seen 
in the SHA can be directly related to the rms multipole 
moments of the peculiar velocity field. 

Equations (9), (10), and (11) are also very similar to 
the harmonics of the temperature fluctuations in the cosmic 
microwave background due to the Sachs- Wolfe effect. 



(la;. 



2c J TT 



dk 
k^ 



PM{k)\ji{2cklH,)f 



(20) 



for Qo = 1 (cL, Bertschinger 1992). Relations 18 and 19 
suggest the possibility of determining, for a given power- 
spectrum in real-space, the product Qo^crg where erg is the 
normalization of the power spectrum of the entire matter 
distribution, independent of biasing. 

In summary, the SHA formalism therefore offers a new 
approach to the study of redshift distortions and the esti- 
mation of Qo which is conceptually rich and allows the unifi- 
cation of a variety of cosmic phenomena, from cosmography 
to quantitative statistical studies, into a single mathematical 
language. 



+ Strictly speaking the bias factor relating the galaxy and mass 
correlation functions is not the same as the bias that enters into 
Equations 1 and 11. Only in the constant bias model assumed 
throughout this paper are the two equal. 



6 Fisher, Scharf, & Lahav 



o.z 



Multiple Window Likelihoods 



— 0.2 

-o 
o 
o 

M-0.4 

CD 



-0.6 



-0.8 



-1 — I — I 1 — I — I 1 — I — I r- 



I 




T" 



r=o.5 



—I I I I L 



T" 



— 1 1 1 1 1 1 1 1 1 1 1 T" 



1.5 



O 
/3 



I 

V 




r=o.2 



_l U I I I I L_ 



1.5 



Figure 3. Likelihood curves (ln£) for f3 for a single mock IRAS realization. The left panel corresponds to likelihood when the correct 
(r = Qoh = 0.5) power spectrum used to construct the mock catalog is used, while the right panel corresponds to the likelihood if the 
model predictions are made with a power spectrum with more large scale power (F = 0.2). The solid curve in each panel is the likelihood 
curve obtained using the value of (Tg derived from the projected harmonics (/(?") =constant) for the assumed power spectrum. The dotted 
and dashed curves correspond to the likelihoods when (Tg is perturbed down and up (respectively) by 10%. The arrows in each panel 
indicate the true value of f3 in the simulation, [3=1. 



50 km s~^. The resulting value of /? was stable to within 
2%. As a further check we divided the simulations into two 
halves and then perturbed the redshifts in one half with a 
Gaussian of dispersion 30 km s~^ and the other half with 
a dispersion of 70 km s~^ in an effort to simulate possible 
systematic errors in the redshifts between the Northern and 
Southern hemispheres; again the derived /? was remarkably 
stable, giving us confidence that the technique is not biased 
significantly by measurement error in the redshifts. 

The correlation function in real space has been deter- 
mined for the 1.2 Jy survey (Fisher et al. 19936). The real 
space correlation function can be integrated analytically to 
derive the variance in a 8 /t~^Mpc sphere erg which is the 
conventional normalization of the power spectrum. Fisher et 
al. (19936) found erg = 0.69 ± 0.04. We have adopted this 
normalization for the power spectrum used in the likelihood. 
The shape of the IRAS power spectrum is somewhat more 
uncertain but appears consistent with a CDM like model 
with F K, 0.2 (Fisher et al. 1993a; Feldmann, Kaiser, & 
Peacock 1993). 

In Table 2, we show the maximum likelihood estimates 
of /? from the 1.2 Jy survey for F = 0.5 and 0.2 with erg 
fixed at 0.69 and when erg is perturbed up and down by its 
1 — (7 errors. If IRAS is well described by a F = 0.2 spectrum 
then the corresponding best value of /? is 0.96; uncertainties 
in the normalization give a 20% variation in this estimate. 
The derived value of /? for a true CDM spectrum (F = 0.5) 



TABLE 2 

IRAS 1.2 Jy i3 estimates 
F = Qo/i 
0.50 0.20 



erg = 0.69 1.34 0.96 
erg = 0.65 1.58 1.16 
erg = 0.73 1.13 0.78 



is higher /? Si 1.3 ± 0.2. The errors in the normalization of 
the power spectrum and the formal uncertainty from the 
likelihood estimator should be independent. Therefore the 
random error in /? is roughly 2^^^0.2 ~ 0.3. The systematic 
uncertainty arising from the shape of the power spectrum is 
< 0.2. 



5 DISCUSSION 

The value of /? derived from the 1.2 Jy sample using the SHA 
technique, /? Si 1.0 ± 0.3, is in agreement with previous esti- 
mates from IRAS redshift surveys. Dekel et al. (1993) found 
/? = 1.281q'59 (2er) by comparing the density field of the 
1.936 Jy IRAS survey (Strauss et al. 19926) with the den- 
sity field inferred from direct measurements of the peculiar 
velocity field using the POTENT algorithm (Bertschinger 

6 Dekel 1989; Dekel, Bertschinger, & Faber 1990). Kaiser 
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a single weighting function, but with a series of different 
weight functions each of which samples a different region of 
space. After experimenting with several different choices of 
weight functions we decided to use four Gaussian windows 
centered at 38, 58, 78, and 98 /t~"'Mpc each with a disper- 
sion of 8 /t"^Mpc. For each of the windows, we computed 
the corresponding weighted redshift harmonics for each of 
the mock IRAS catalogs. Next, for a given choice of power 
spectrum, we maximized the likelihood given in Equation f 4 
over a grid in /?. Since our analysis is valid only in the linear 
regime, we restricted the likelihood computation to I < fO. 

The likelihood of measuring a set of redshift harmonics 
with N weighting functions, {i = 1, . . . N] is given by 

l„/:^_(?l±l)l„|A|_l^A-i [Reial,)Reiai„)] (f4) 

i,j m=l 

where Atj is the covariance matrix 

A^J = {aima\*^)TH + {aima\*^)sN , (f5) 
with 

CO 



(^,«o)(fc)+^^,-o)(fc)y (16) 

and 

CO 

{al^aZ)sN = j drr'Hr)[/'\r)/^^\r)] , (f7) 



corresponding to the different weighting functions, f^'^s). 
In Equation f4 \A\ is the absolute value of the determinant 
of A,j, and Re(aim) and Im(aim) refer to the real and imag- 
inary parts of aim (SL93). The superscripts on ^'"(A;) and 
^'^(A;) refer to the window function of the corresponding 
weighting function. Since the harmonics at different I are 
independent, the total likelihood for I < Imax is just given 
by the sum of the likelihoods in Equation f4 for each I up 
to Imax- For our choice of weight functions, the off-diagonal 
elements of the covariance matrix, Atj, are ^ fO% of the 
largest diagonal elements; the four spatial windows defined 
by the weighting functions therefore give roughly indepen- 
dent estimates of /?. 

Table f shows the resulting estimates of /? for the nine 
mock catalogs when the correct normalization and shape of 
the power spectrum were used in the likelihood. The mean 
value of /? for the nine mock samples is 0.94 ± 0.2 which is 
in excellent agreement with the actual value /? = f. This 
quantifies the qualitative agreement between the predicted 
and observed distortion seen in Figure 2; when the correct 
power spectrum is used the likelihood estimate for /? is un- 
biased. The formal uncertainty in a single realization (given 
by A£ = — f/2) is comparable to the scatter over the nine 
realizations, i.e.. A/? = 0.2. 

How do the estimates of /? change when the incorrect 
power spectrum is used in the likelihood? Figure 3 shows the 



TABLE 1 

CDM MONTE CARLO RESULTS 



Realization /? 

f 0.93 

2 f.ff 

3 f.06 

4 0.86 

5 0.83 

6 0.88 

7 f.fS 

8 f.f3 

9 0.5f 
Mean: 0.94 



Standard Deviation: 0.20 



likelihood as a function of /? for a single mock IRAS cata- 
log when the shape and normalization of the power spec- 
trum are altered from their true values. The left hand panel 
shows the result of varying the normalization of the true 
power spectrum by ±fO%. When the normalization is too 
low by fO% the corresponding maximum likelhood value of 
/? is systematically too high by roughly 40%; the converse 
is true when too high a normalization is used. In the right 
hand panel of Figure 3 we show the values of /? for a power 
spectrum with more large scale power (F = 0.2) than used 
in the simulations. This power spectrum predicts excessive 
power on large scales and therefore the best fit value of /? is 
systematically low. 



4 Qo FROM THE 1.2 JY IRAS SURVEY 

We now proceed to apply the formalism described in the pre- 
vious section to redshift survey of 5313 IRAS galaxies flux- 
limited to 1.2 Jy at 60/tm, selected from the IRAS database 
(Strauss et al. 1990; Fisher 1992). The sample covers 87.6% 
of the sky and is complete for \h\ > 5° with the exception 
of a small area of the sky not surveyed by IRAS . The SHA 
analysis discussed above relies on complete Itt steradian sky 
coverage. Although statistical corrections can be applied to 
the harmonics with partial sky coverage (Scharf et al. 1992; 
SL93), we have adopted a simplier method of dealing with 
incomplete sky coverage. We have interpolated the redshift 
data through the plane in a way which smoothly contin- 
ues structure; this is the same approach used by Yahil et 
al. (1991) in their reconstruction of the IRAS peculiar ve- 
locity field. Moreover, regularized mask inversion of the har- 
monics shows the effect of incomplete sky coverages is small 
for the IRAS 1.2 Jy sample geometry and I < 10 (Lahav 
et al. 1993, in preparation). The selection function, 4>{r), 
is computed using the techniques of Yahil et al. (1991). 
In practice the selection function is computed by assigning 
luminosities based on galaxy redshifts, not distances; fortu- 
nately, the shape of the 4>{r) is insensitive to the flow model 
used to correct for the effects of peculiar velocities (Saunders 
et al. 1990). 

The typical redshift error in the 1.2 Jy IRAS survey 
is ~ 50 km . We repeated the TV-body tests discussed 
in the previous section with the redshifts of the simulation 
particles perturbed by a Gaussian deviate with dispersion 
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(1987), Kaiser & Lahav (1988), and Strauss et al. (1992a). 
We stress, however, that in linear theory the choice of ref- 
erence frame (i.e., the value of Vobs) affects only the value 
of the dipole harmonic and that for I > 2 the results are in- 
dependent of the assumed motion of the observer. The rms 
averaging procedure in Equation (11) treats the observer as 
a typical point in space and therefore assumes that the ob- 
server's peculiar motion is not unusually high or low. This 
point should be kept in mind if the motion of the Local 
Group is atypical. In practice, these complications can be 
avoided by working in the comoving frame of the observer, 
i.e., the reference frame in which Vobs vanishes. 



Mock IRAS N-body Test 



TV-BODY TESTS: 



A WAY TO DETERMINE 



In order to check the validity of our formalism, we computed 
both the real and redshift space weighted harmonics in an 
TV-body simulation of a standard Cold Dark Matter uni- 
verse characterized by QqA = 0.5. The simulations evolved 
64'^ particles in a box of 180 /t~^Mpc using the P"^ M algo- 
rithm until rms variance of the density field in a sphere of 
8 /t~^Mpc reached ag = 0.61; further details of the simu- 
lations can be found in Frenk et al. (1990). In an effort to 
mimic current observational data, we extracted nine mock 
galaxy catalogues designed to closely match the properties 
of the 1.2 Jy IRAS redshift survey (cf., Strauss et al. 19926, 
Fisher 1992). Galaxy candidates were chosen as unbiased 
tracers of the underlying particle distribution and therefore 
/? = 1 in the mock catalogues. The procedure for extract- 
ing these catalogues is described in Gofski et al. (1989) and 
Fisher et al. (1993a). 

Figure 2 shows the mean weighted harmonic power 
spectrum, Cf , defined as the azimuthally averaged harmonic 
by the shot noise contribution to the harmonic (cf., Scharf 
et al. 1992), 



21 + 1 



+ 1 

E 



a;^ 



cf 



(13) 



The harmonics were computed both in real and redshift 
space by performing the direct summation given in Equa- 
tion (4), squaring the resulting harmonic and then subtract- 
ing the shot noise contribution as given in Equation (12). 

In deriving the harmonics shown in Figure 2, we used a 
somewhat ad hoc Gaussian weighting function for f{s) with 
centroid at 40 /t~^Mpc and dispersion a = 15 /t~^Mpc. The 
redshifts in the simulation were corrected to the comoving 
frame of the observer, i.e., the frame where "Vobs • r = 0. The 
solid symbols in Figure 2 show the mean value of the real 
space harmonic spectrum computed from the harmonics in 
each of the nine mock IRAS catalogues. The open symbols 
correspond to the same calculation performed in redshift 
space. 

The error bars correspond to the variance in the har- 
monic estimates over the nine simulations. The solid and 
dashed curves correspond to the theoretical prediction given 
by Equation (11), using the known real space power spec- 
trum and P of the simulations. From Figure 2, we see that 
for our choice of window function, the redshift space har- 




Figure 2. The harmonic power spectrum in real space (filled 
symbols) and redshift space (open symbols) for the mean of nine 
mock IRAS catalogues. The solid and dashed lines represent the 
linear theory predictions (cf., Equation (11)) for real and redshift 
space respectively. 



monies are enhanced and that the overall amplification is 
described quite well by the prediction of Equation (11). 

Since the distortion is well described by our formalism, 
the next step is to ask whether one can predict the value 
of P given a set of harmonics measured from a redshift cat- 
alogue. The predicted distortion depends on knowing the 
correct shape and amplitude of the galaxy power spectrum. 
The normalization of the power spectrum can be derived by 
looking at the amplitude of the number weighted (/(s) = 1) 
harmonics since in that case the redshift distortion is iden- 
tically zero; this was the approach adopted by SL93. Alter- 
natively, the amplitude of the power spectrum in real space 
can be fixed by looking at the real space correlation function. 
As we will see in the next section, this approach is feasible 
for redshift surveys large enough to allow the real space cor- 
relation function to be deprojected from the redshift space 
correlation. 

The shape of the power spectrum can be conveniently 
parametrized by a series of phenomological CDM models 
with with varying F = QqA (e.g., Efstathiou, Bond, & White 
1992). Of course, in the TV-body simulations we know the 
power spectrum precisely; however in order to apply the 
formalism to actual data we must quantify the dependence 
of the derived /? on the shape of the power spectrum. 

The weighting function, f{r), used to construct the har- 
monics can be chosen arbitrarily. Clearly, one would like to 
pick the weight function which simultaneously maximizes 
the distortion while minimizing the noise. We are currently 
working on this optimization problem (cf., § 5 below). How- 
ever, the arbitrariness of the weighting function can be some- 
what circumvented by computing the harmonics not with 
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Figure 1. An example of the window functions appearing Equations (9) and (10). In each panel, the dotted curves show k'^P{k) for 
the CDM {p,oh = 0.5) linear theory power spectrum and the dot-dashed curves show the corresponding CDM nonlinear power spectrum 
determined from the A^-body simulation described in the text. In the left panel, the solid curve represents the |^^(A;)p, the window 
function corresponding the real space contribution to the quadrupole (/ = 2) moment as defined in Equation (9). The solid curve in 
the right hand panel respresents the window function of the redshift distortion (again for / = 2), [^^(fc)^ defined in Equation (10). 
In deriving the window functions we took the weighting function, to be a Gaussian with mean of 40 /i~"^Mpc and a dispersion of 

15 /i~"^Mpc; this corresponds to the window used construct the harmonics in Figure 2. 



In Equation (7) and below, j'i(kr) = dji(kr) /d(kr) refers to 
the first derivative of the Bessel function. In Equation (7), 
we have used the continuity equation to relate the velocity 
to the density field, Vj, = — i/95j!k/A;^ ; in this limit Equa- 
tion (6) can be written (for I >1), 

afm = ^/ d'^^k [*f(fc)+/?*r(fc)] l^„(k) , (8) 
where 

CO 

*f(A;) = j drr^cj>{r)f{r)ji{kr) (9) 



describes the real space contribution to the harmonics and 

CO 

^f{k) = \j dTT'4>{r)^ (][{kT) - , (10) 



is a "correction" term which embodies the redshift distor- 
tions. In Equation (10), Snisa, Kronecker delta which con- 
tributes only to the dipole (I = 1) harmonic. In the case of 
equal weight, df / dr = 0, '$f(k) vanishes and Equation (8) 
leads to the real space expression for the harmonics derived 
by Scharf et al. (1992). The expressions given in Equa- 
tions (8), (9), and (10) are exact in linear theory. 

The window functions, ^'f (A;) and '^f(k), correspond- 
ing to a Gaussian weighting function /(s) with mean 40 
/t~^Mpc and dispersion 15 /t~^Mpc are shown in Figure 1. 
The distortion term, '^f (k), (the right hand panel) is sensi- 
tive to the nature of the velocity field and therefore to the 



long wavelength modes of the density field; in contrast, the 
real space window, ^'f (A;), is relatively insensitive to large 
scale fluctuations. 

The expected linear theory rms value of the harmon- 
ics is given by ensemble average of the square of Equa- 
tion (8). Using the statistical independence of Fourier 

wavest, {S^S^) = (27r)^P„(A;)«(^)(k - k'), this yields 

CO 

{\aLf)=l J dkk'Pn{k)\-fr{k) + l^-f?{k)\' . (11) 



In real data the square of harmonics in Equation (11) will 
have a discreteness or "shot" noise contribution; this can be 
modeled by adding 

CO 

{\ai^f}sM = J rfrr^0(r)[/(r)f (12) 



to Equation (11). 

The I = 1 or dipole distortion has two contributions: one 
from the external dipole moment of the velocity field around 
the observer and one induced by the motion of observer it- 
self. The latter distortion is caused by "Vobsi^ and is the 
origin of the so-called "rocket effect" discussed by Kaiser 

^ This assumes that the survey contains many independent 
modes of the wavenumber in question and usually phrased some- 
what loosely as the "fair" sample hypothesis. 
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In practice there are several disadvantages to analysing 
the redshift distortions with {s(s). Firstly, Kaiser's analy- 
sis is valid only in the linear regime where 8{r) ^ 1. Un- 
fortunately, {s(s) depends on cumulative moments of the 
real space correlation function which can contain contribu- 
tions from nonlinear evolution on small scales (although see 
Hamilton et al. 1991 for a discussion of moments which are 
unaffected by nonlinear evolution). Secondly, Equation (1) 
is strictly only applicable to volume limited samples. In the 
case of a flux limited sample characterized by a radial selec- 
tion function, 4>{r), (the number density in a homogeneous 
flux limited survey being oc <j)[r)) there is a contribution 
to the redshift distortion which is given by 



2 + 



U{r) - Vois 



(3) 



where we have adopted Kaiser's notation for the radial com- 
ponent of the peculiar velocity field, U{r) = v(r) -r and Vobs 
is the velocity of observer defined in the same reference frame 
as v(r). Although Kaiser considered the effect of a radial se- 
lection function in his analysis, this term was assumed to be 
small and was omitted from Equation (1). 

In this paper, we present a new and complementary way 
of investigating redshift distortions based on a linear ex- 
pansion of the galaxian density field in spherical harmonics. 
Spherical Harmonic analysis (hereafter SHA) leads to a sur- 
prisingly simple expression for the redshift distortion which 
naturally includes the effect of a radially varying selection 
function. Moreover, since the l"^ order harmonic probes a 
typical scale irD/I (where D is roughly the sample depth), 
one can ensure the analysis is carried out in the linear regime 
by restricting the investigation to the low order multipole 
moments. 

Scharf et al. (1992), Lahav et al. (1993), and Scharf & 
Lahav (1993; hereafter SL93) have shown that SHA pro- 
vides a useful method for probing the local cosmography of 
the galaxy distribution. During their investigation of galaxy 
clustering using SHA, they noted the sensitivity of their re- 
sults to assumptions regarding the nature of the galaxy pe- 
culiar velocity field. The present work is an extension of the 
SHA which statistically incorporates the effect of redshift 
distortion in a self-consistent way. 

We begin the paper by outlining the SHA method in § 2. 
In § 3, we check our formalism by comparing our predictions 
for the harmonics in redshift space with those actually com- 
puted from numerical TV-body simulations and show how it 
can be used to formulate a maximum likelihood estimator 
for /?. In § 4, we apply our formalism to the IRAS 1.2 Jy 
redshift survey and derive an estimate of /?. We conclude in 
§ 5. 



2 METHOD 

Following SL93, we define a weighted spherical harmonic 
decomposition of the flux-limited density field, in redshift 
space as. 
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where Ng is the number of galaxies in the survey, Yim is the 
usual spherical harmonic, and f{s) is an arbitrary radial 
weighting function. In the special case of constant weight- 
ing all redshift distortions vanish and one simply obtains 
the projected angular distribution on the sky. The harmon- 
ics derived from such weighting can, for example, be used to 
derive the real space correlation function, {(r), using tech- 
niques developed by Peebles (1973). However, in the more 
general case where f{s) varies with s, peculiar velocities will 
distort the harmonics. 

In order to relate the redshift harmonics given by Equa- 
tion (4) to their real space counterparts, we first rewrite the 
summation in Equation (4) as a continuous integral over the 
density fluctuations in redshift space. 
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(5) 



In Equation (5), the selection function is assumed to be nor- 
malized such that J dr 4'{r) = Ng/uj where uj is the solid 
angle subtended by the survey. Notice that the selection 
function in Equation (5) is evaluated at the galaxy's dis- 
tance, not redshift, because if the catalogue is flux limited 
the probability of a galaxy being at redshift, s, will be pro- 
portional to the selection function evaluated at the galaxy's 
actual (albeit unknown) distance, i.e. oc i/i(r). Next, we note 
that, by construction, ras(s)rf'^s = nn^rjd'^r, where ras(s) 
and raR(r) refer to the densities in redshift and real space 
and that if the perturbations induced by peculiar motions 
are small, then we can perform a Taylor series expansion of 
all redshift quantities to first order in the density fluctuation, 
e.g., /(s) ~ /(r) -I- {U{r} - Yobs ■ r). Thus, expansion 
of Equation (5) yields. 
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d'r(t>(r)f(r)[l + Snir)]YUr) 
d'r <t,{r}^ {U{r} - Yobs ■ r}Yi,^{f} 



(6) 



There are several comments to be made regarding the 
validity of Equation (6). Firstly, we have assumed that the 
summation in Equation (4) is carried out over all galaxies 
in a flux limited redshift survey. In this case, the integrals in 
Equation (6) extend over all space and there are no "surface" 
terms arising from the deformation of the boundary of the 
integration region that occurs in the transformation from 
redshift to real space. Secondly, there is an apparent absence 
of terms involving the derivatives of the selection function; 
these terms would be manifest if we, like Kaiser (1987), took 
f{s) to be the special case of l/i/i(s). 

The expansion in Equation (6) can be simplified by ex- 
panding the radial peculiar velocity field dependence, U{r), 
in spherical harmonics: 



C/(r) 



v(r) • r 



'^jUkr)YiUmU^) ■ (7) 



(4) 



* In our analysis f{s) is required to be continuous in its first 
derivative and to vanish at infinity; this simplifies the analysis by 
eliminating surface terms that arise when f{s) has a discontinuous 
boundary. 
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ABSTRACT 

We examine the nature of galaxy clustering in redshift space using a method based 
on an expansion of the galaxian density field in Spherical Harmonics and linear theory. 
Our approach provides a compact and self-consistent expression for the distortion when 
applied to fiux limited redshift surveys. The amplitude of the distortion is controlled 
by the combination of the density and bias parameters, j3 = Q'^ '^/h; we exploit this fact 
to derive a maximum likelihood estimator for /3. We check our formalism using #-body 
simulations and demonstrate it provides an unbiased estimate of /3 when the amplitude 
and shape of the galaxy power spectrum is known. Application of the technique to the 
1.2 Jy IRAS redshift survey yields /3 1.0; both random errors (from counting statis- 
tics and the uncertainties in the power spectrum normalization) and systematic errors 
(from the uncertainty in the shape of the power spectrum) individually contribute 20% 
uncertainties in this estimate. This estimate of /3 is comparable (both in amplitude and 
uncertainty) with previous measurements based on comparisons of the IRAS density 
field with direct measurements of peculiar velocities and analyses of the acceleration of 
the Local Group, but the Spherical Harmonic Analysis has the advantage of being easy 
to implement and is largely free of systematic errors. 

Key words: Cosmology: theory-large-scale structure 



1 INTRODUCTION 

The clustering of galaxies in redshift space appears system- 
atically different from the clustering that one would observe 
in real space. As Kaiser (1987) has pointed out, structures 
tangential to an observer's line-of-sight will appear enhanced 
in redshift space due to the coherent motions of galaxies on 
large scales. Using linear theory. Kaiser was able to deduce a 
simple expression for the resulting direction dependent am- 
plification of the redshift space density field. Kaiser found 
that the redshift space power spectrum was anisotropic and 
related to the real space power spectrum (assumed to be 
isotropic) by 

Ps(]^) = Pn(k) (l + /?/tL)' , (1) 

where (Ikl is the cosine of the angle between the wave vector 
k and the observer's line-of-sight. In Equation (1) and in the 
remainder of this paper, we adopt the notation /? = 
where Qo and b are the current density and bias (assumed 
here and below to be independent of scale) parameters re- 
spectively; we also will use R and S super and subscripts 
when referring to quantities in real and redshift space re- 
spectively. The quantity /? naturally appears in expressions 



for the redshift distortion since it is the proportionality fac- 
tor between the velocity and density field in linear theory, 
V • V = -HoiSS (e.g., Peebles 1980). 

While the redshift distortion given by Equation (1) is 
quite simple in the Fourier domain, it leads to a fairly com- 
plicated expression when cast in terms of the redshift space 
correlation function, {s(s). The full angular dependence of 
{s(s) was first derived by Lilje & Efstathiou (1989) and then 
by McGill (1990) and Hamilton (1992) in various forms. We 
find it illuminating to write {s(s) in the form 

is(s) = [l + + 1/?^] Ao(s)To(jiLs) 

- 11^+^^'] Ms)Mt^Ls)+ [^/?'] A,{s)v,i^l,s) , (2) 

where Ai(s) = J dk Pn^kjji^ks), ji{x} is the spherical 
Bessel function and Vigils) is the usual Legendre polyno- 
mial evaluated at the cosine of the angle between s and the 
observer's line-of-sight. In principle, the dependence of the 
redshift distortions on /? offers the hope that Qo can be de- 
termined dynamically (e.g., Hamilton 1992, 1993) assuming 
the bias can be determined independently. 



